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ABSTRACT 



Aims. We aim to exploit the available INTEGRAL/SPl data to provide time-averaged spectra of the brightest soft gamma-ray sources. 
Methods. Employing a maximum-likelihood fit technique for our SPI data analysis, we take as input to our source model the source 
catalog derived by Bouchet et al. (2008) from a SPI all-sky study. We use the first four years of public SPI data and extract spectra 
between 25 keV and 1 MeV for the 20 catalog sources detected by Bouchet et al. at 200 - 600 keV with > 2.5 tr. In order to verify 
our analysis, we also extract spectra for the same sources from the corresponding INTEGRAL/lSGRl data. We fit adequate spectral 
models to the energy range 25-1000 keV for SPI and 25-600 keV for ISGRI. We use our spectra from the Crab (which is among the 
20 sources studied here) and an empty location in a crowded field to derive an estimation of the systematic errors. 
Results. The agreement between our SPI and ISGRI measurements is good if we normalise them on the Crab spectrum. Our SPI flux 
measurements also agree well with those by Bouchet et al. (2008). All 20 sources in our sample are detected independently in the 
bands 25-100 keV and 100-200 keV. At 200-600 keV we detect eight sources, at 600-1000 keV we detect two sources. Our spectra 
agree well with the results from previous publications where available. For six of the 14 XRBs in our sample we find evidence for 
a hard powerlaw-component which becomes dominant above the cutoff energy of the thermal Comptonization component. In two of 
these cases, our study provides the first indication of such emission. For the others, our results confirm previous studies. Our spectrum 
of the Crab, integrated over 1.3 Ms, shows a significant flux in all points and is well described by a powerlaw with a break near 
100 keV and spectral indices 2.1 1 and 2.20. 

Key words. Gamma-rays: observations - X-rays: individuals: Crab, Vela Pulsar, NGC 4151, NGC 4945, Cen A, XTE J1550-564, 
4U 1630-47, Swift J1656.3-3302, OAO 1657-415, GX 339-4, 4U 1700-377, IGR J17091-3624, GX 354-0, IE 1740.7-2942, IGR 
JI7464-32I3, GRS 1758-258, Ginga 1826-24, GRS I9I5-HI05, Cyg X-I, Cyg X-3 - X-rays: binaries 



1. Introduction 

Soft gamma-rays have to be observed with satellite-based in- 
struments, and so far there have been only three major mis- 
sions taking data in the range from 100 keV to a few MeV 
with adequate sensitivity: GRANATwith the SIGMA instrument 
(Revnivtsev et al. 2004), CGRO wit h OSS E (Johnson et al. lT993T l 
and COMPTEL (Sc honfel der et al. l2000] l. and INT EGRAL with 
IBIS (Ubertini et al. 12003] ). and SPI (Vedrenne et al. 12003] ). 

On INTEGRAL, the only soft gamma-ray observatory oper- 
ational today, both the high-energy instruments, IBIS and SPI, 
are capable of measuring point source spectra. They comple- 
ment each other in that IBIS provides angular resolution and 
good sensitivity at energies below a few 100 keV while SPI has 
superior energy resolution and high-energy sensitivity and per- 
mits the study of extended diffuse emission with its large field of 
view. 

Sources of soft gamma radiation with photon energies be- 
tween 100 keV and several MeV are astrophysically interest- 
ing because they permit us to study matter at supra-thermal en- 
ergies. Essentially all known point-like emitters of such radia- 
tion consist of compact objects such as black holes or neutron 
stars which convert the gravitational energy of their surrounding 
matter or their own extreme rotational and magnetic energy into 
translational kinetic energy and then in turn partially into high- 
energy photons. The black holes come either as stellar black 



holes in X-ray binaries or as supermassive black holes in the cen- 
tres of galaxies. All other sources of soft gamma-rays are spa- 
tially extended or unresolved and high in number-density. They 
form a continuum which seems, in the case of our Galaxy, to ex- 
tend over the bulge and a significant part of the disk (Bouchet et 
al. 12005] Revnivtsev et al. .2006. Krivonos et al. 2007 ). 

In the most recent all-sky st udy o f INTEGRAL/SPl soft 
gamma-ray data, Bouchet et al. (I2008I I have presented a de- 
tailed analysis of both diffuse and point source emission in the 
energy range from 25 keV to 600 keV based on all SPI data 
from February 2003 up to May 2006. They conclude that the 
diffuse component becomes more prominent with increasing en- 
ergy reaching 32 % in the central radian above 100 keV. The 
diffuse emission is bright over up to +45° in galactic longitude 
and ±10° in galactic latitude depending on energy. 

Bouchet et al. (2008 ), BOS in the following, also present sep- 
arate catalogs of point sources for the four energy ranges 25-50 
keV, 50-100 keV, 100-200 keV, and 200-600 keV which they 
detect at least at a significance level > 2.5cr. As the sensitivity 
of SPI relative to the steeply dropping source spectra decreases 
with increasing energy and the diffuse emission becomes more 
dominant, the number of detected point sources decreases with 
increasing energy as well. In the energy range 200-600 keV a 
total of 20 sources are detected. 

In this work, we take the BOS catalog for the energy range 
200-600 keV and derive time-averaged INTEGRAL/SPl spectra 
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for each of the 20 point sources contained in it. The resuhing 
catalog of spectra of the 20 brightest soft gamma-ray sources 
will be useful to, among others, modellers of the hard X-ray and 
soft gamma-ray background caused by unresolved sources of the 
same source classes as in our catalog, high-energy observers of 
these sources who have to average over longer time-scales in or- 
der to achieve detections with acceptable significance, and gen- 
erally researchers interested in the energy balance of the galaxy 
(most of the 20 sources are galactic). 

We employ a newly developed maximum likelihood fit anal- 
ysis technique based on the spimodfit software by Halloin & 
Strong ^2MT\i to all public INTEGRAL data available for these 
sources (i.e. up to November 2006, the status when our dataset 
was frozen). In order to minimise inhomogeneities of our sys- 
tematic uncertainties, the analysis method is not optimised for 
any individual source, in particular not for the Crab. 

For comparison and verification of our method, we also de- 
rive the spectra for the same sources from the corresponding 
INTEGRALIimSllSGm data. 

This paper is meant to serve two purposes: (a) be a supple- 
ment to BOS and investigate the average spectra of the brightest 
known soft gamma-ray sources, and (b) demonstrate an analysis 
technique for INTEGRAL/SPl data based on spimodfit. 

In Sect. |2] the SPI and ISGRI instruments are briefly intro- 
duced, and the observations and the data selection are described. 
Section [3] describes the analysis technique. Section |4] presents 
and discusses the spectra individually while Sect. |5] summarises 
the results, discusses them jointly and concludes. 

2. Instruments and observations 

This study uses data from the SPI and the ISGRI instruments on 
INTEGRAL. The properties of these instruments are discussed 
in detail elsewhere (see Vedrenne et al. (I2003I I and Roques et al. 
( l2003T l for SPI and Ubertini et al. ([2003T l for ISGRI). Here we 
just briefly summarise the main aspects. 

The ISGRI instrument is part of the high-sensitivity coded- 
mask imager IBIS. It has a fully coded field of view of 9°x9°, 
an angular resolution of 12' FWHM, and a spectral resolution of 
9% at 100 keV. Its well calibrated energy range is 14 keV to ca. 
700 keV. The nominal continuum sensitivity (AE = E/2, 3 cr in 
10^ s) at 100 keV is 3 mCrab. 

The SPI instrument is the coded-mask imaging high- 
resolution spectrometer on INTEGRAL. Its fully coded field of 
view has a radius of 8°. With its 19 pixels (hexagonal, cooled 
Ge detectors) it achieves an angular resolution of 2.5° FWHM. 
The average spectral resolution of the detectors deteriorates only 
marginally as a function of energy from AE a; 1.8 keV below 
200 keV to Afi a; 4 keV at 3 MeV. This means, the relative en- 
ergy resolution improves from 0.9 % at 200 ke V to 0. 1 3 % above 
3 MeV. The calibrated energy range is 25 keV - 8000 keV. The 
nominal continuum sensitivity (as defined above for IBIS) at 100 
keV is 16 mCrab. But above 1 MeV, the continuum sensitivity 
for point sources becomes inadequate. Only the Crab is detected 
at 1-8 MeV, and an electronic noise problem between 1 and 2 
MeV results in large systematic errors in this interesting energy 
range. We therefore limit this study to the range 25-1000 keV. 

Both energy resolution and sensitivity are time dependent; 
Radiation damage degenerates the detectors leading to a deteri- 
oration of the energy resolution by a few tenths of a keV (neg- 
ligible for continuum analysis) over timescales of months. This 
degradation is cured by regular annealings (controlled heatings) 
of the detectors. An annealing period lasts typically 12 days dur- 
ing which no science data can be taken. 



While the spectral resolution is hence kept constant for our 
purposes, the sensitivity of SPI remains time dependent because 
the continued exposure to cosmic radiation leads to increasing 
background from the activation of all parts of the INTEGRAL 
satellite. In addition, the intensity of the cosmic radiation itself 
is strongly time dependent because of solar wind variablity (so- 
lar cycle) and solar flares. Over the course of the INTEGRAL 
mission so far, this has lead to a nearly steady increase in instru- 
mental background. In November 2006, the instrumental back- 
ground had increased by about 80% compared to the time shortly 
after launch (November 2002). It is expected that this trend will 
not continue (and possibly be reversed) beyond the end of the 
present solar minimum which has been a very prolonged one. 

2.1. Data selection 

The INTEGRAL data, as they are delivered to scientists, are sub- 
divided into time periods of constant pointing, so-called science 
windows which typically have a duration of about 30 minutes. 
A second relevant time unit is the "revolution" or orbit which 
is important to the science analysis because the highly excen- 
tric orbit of the INTEGRAL spacecraft takes the instruments 
through the Earth's radiation belts every three days ("perigee 
passage") making it necessary to temporarily switch off many 
of the INTEGRAL instruments including ISGRI. Observation 
scheduling and public data release is then organised per revo- 
lution. 

This study is based on all available high-quality public data 
(at the time the dataset was frozen by us). These are the data 
from revolutions 21 to 500, i.e. from 15 December 2002 to 18 
November 2006. From this time period, we select good data for 
each of the objects for which we want to derive spectra. The 
selection criteria are the following: 

1 . The angular distance between the pointing direction and the 
position of the object of interest is required to be less than 
10°. This is a compromise between maximum exposure time 
and maximum data quality. For SPI the maximum off-axis 
angle could be further increased, but for ISGRI, accepting 
data beyond 10° radius (which is well into ISGRJ's partially 
coded field of view) reduces the accuracy of the spectra sig- 
nificantly. Since INTEGRAL does mainly pointed observa- 
tions based on individual proposals (which have emphasised 
galactic objects so far), the sky exposure achieved by the 
INTEGRAL instruments up to now is far from uniform. At 
low galactic latitudes our dataset therefore contains up to a 
few thousand science windows while only a few hundred are 
available for the typical object at higher latitudes. 

In those cases where a very large exposure is available, we 
can reduce the selection radius to 8° in order to improve 
data quality. Furthermore, if several objects of interest form a 
closely spaced group, we use a single dataset centred roughly 
on the centre of the group. 

2. Time periods where the data quality was known to be 
bad such as times near the radiation belt entry or exit of 
the spacecraft, times shortly after annealings, times during 
which solar flares strongly raised the background, or times 
during which hardware problems occurred, were excluded 
based on the publicly available database of such events in- 
cluded in the OSA software (version 7.0) provided by the 
INTEGRAL Science Data Centre (ISDC, Courvoisier et al. 
120031) . This results in the removal of about 31 % of the 
archival data (based on elapsed time and averaged over the 
entire mission up to revolution 500). 
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3. In order to exclude time periods with problems which might 
have been missed in the compilation of the database used in 
the previous criterion, quality cuts were imposed on a num- 
ber of "housekeeping" parameters using the SPI housekeep- 
ing database compiled at MPlfl In particular, the cuts lim- 
ited the values of the SPI detector temperature, the veto rate 
of the SPI Anti-Coincidence Shield (ACS), the rate of SPI 
single interaction events, and the orbital phase. The cuts were 
derived from inspecting the long term behaviour of these pa- 
rameters and assuming smoothness of their time dependence. 
The orbital phase was conservatively limited to the range be- 
tween 0.1 and 0.87 in order to safely exclude any radiation 
belt influence. These housekeeping parameter cuts removed 
another 9.5 % of the archival data (percentage defined as 
above), mostly at the beginning of the mission. 

For detailed analysis of the 20 objects from the BOS catalog 
we finally utilise the datasets as presented in Table [T] 

3. Data analysis 

3.1. SPI data analysis using spimodfit 

The spectral data analysis was carried out using the software 
packages spiselectscw v3.9 and spimodfit v3.0 (both written by 
Halloin & Strong 2007) with additional infrastructure software 
written by DP. The spimodfit package is part of the standard 
INTEGRAL data analysis software (OSA) since version 7.0 . 

Using spiselectscw, the data selection criteria described in 
Sect. l2.1l were applied to the calibrated and pre-binned data from 
the SPI database at the Max Planck Institute for extraterrestrial 
Physics. Using spimodfit, a source region model was then fit- 
ted to the selected data. The parameters of this model are in our 
case dynamically constructed by spimodfit from two main com- 
ponents: 

1 . an isotropic background component. 

2. the contributions from those point sources which are within 
20° of the optical axis of the instrument for a given science 
window. 

The isotropic component was modelled as the product of a com- 
mon scaling factor and a relative sensitivity for each of the 19 
detectors of SPI (taking into account the two detector failures 
during the course of the mission so far). The relative detector 
sensitivity was permitted to vary on a timescale of 10 revolutions 
(nearly 30 days) as a compromise between assuming no variabil- 
ity at all (resulting in a degraded fit) and the minimum reasonable 
variability timescale of one revolution (approx. 3 days, resulting 
in larger uncertainties due to the increased number of fit param- 
eters). It was verified that assuming a shorter timescale than 30 
days did only change the results within their statistical errors. 

The common background scaling factor was permitted to 
vary on a short timescale of 2 h (corresponding to typically four 
science windows). This was the conservative choice determined 
by the actually observed changes in the overall event rate (which 
is strongly dominated by background) and the aim to minimise 
the number of model parameters in order to improve sensitivity. 

The model components stemming from the point sources are 
redetermined for every pointing: The algorithm first determines 
which sources from the input catalog are within 20° of the point- 
ing axis. Then the tabulated detector response (taking into ac- 
count the absolute time, spacecraft orientation, the coded mask 

' |http://www.mpe.mpg.de/gaimna/mstrumeiits/mtegral/spi/w\vw/public-ciata| 



pattern, and the general spacecraft and instrument mass model) 
is used to calculate the relative contribution of each source to 
the counts in each detector. The model fit parameter is the com- 
mon scaling factor of the contributions in each detector, i.e. the 
source flux in the time window corresponding to the pointing. 
Depending on the variability timescale assumed for the individ- 
ual source, this flux is then assumed constant for a certain num- 
ber of pointings. 

When the model is constructed, all parameters are tested for 
relevance, and parameters which do not show an effect on the fit 
function are excluded. Using a maximum likelihood fit which is 
initialised from a ;if^-fit, the model is then fit to the data and the 
statistical errors are determined from the covariance matrix. For 
more details see Halloin & Strong (|2007| l. 

The entire analysis is carried out independently in each 
energy bin such that boundary conditions such as variability 
timescales and the input point source catalog can in principle 
be chosen diff'erently for each energy bin. 

In the largest of our datasets (the one containing 
GRS 1758H-258, see Table[T]l, the number of degrees of freedom 
is 74292 per energy bin, while the number of fit parameters is 
2582. In the smallest dataset (the one containing NGC 4151) 
there are only 665 degrees of freedom and 38 fit parameters per 
energy bin. Thus, for an entire SPI spectrum with 17 energy bins 
(see also the next section), the total number of degrees of free- 
dom in our datasets is between « 11000 and a; 1.3 million, the 
number of fit parameters between 600 and 44000. 

The number of fit parameters is dominated by the back- 
ground scaling factors and relative detector sensitivities. The 
source flux for an individual source from the input catalog con- 
stitutes only one fit parameter for each energy bin. 

Unless otherwise noted, the reduced values of the model 
fit in all energy bins for all results presented here were very close 
to 1.0. 

3.1 .1 . Energy binning 

Our default energy binning was chosen a priori based on consid- 
erations of SPI energy resolution and flux sensitivity. It consists 
of 16 energy bins of logarithmically increasing width between 
25 and 1000 keV. In addition, a narrow energy bin around 511 
keV was inserted in order to check for possible annihilation ra- 
diation features. 

As the analysis showed, many of our sources of interest are 
too weak to be detected significantly above 200 keV in individ- 
ual bins of even this coarse (compared to the SPI energy reso- 
lution) energy binning. So a second wider binning consisting of 
1 1 bins was introduced for weaker sources. The narrow and the 
wide binning have the first three bins (up to 48 keV) in common, 
then the wider binning roughly combines every two narrow bins 
up to 500 keV and after that every 3 narrow bins. The bin bound- 
aries of the narrow binning are 25, 31, 39, 48, 60, 75, 93, 116, 
144, 179, 223, 278, 346, 502, 520, 668, 832, and 1000 keV. The 
bin boundaries of the wider binning are 25, 31, 39, 48, 70, 103, 
150, 219, 320, 502, 520, and 1000 keV. 

3.1.2. Input catalog construction 

As mentioned above, we base our input catalog on the one com- 
piled by B08. This catalog comprises in total 173 sources and 
states fluxes for them in four energy bands (25-50 keV, 50- 
100 keV, 100-200 keV, and 200-600 keV). At low energies, all 
sources of the Bouchet catalog have to be included in the in- 
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Table 1. The 20 objects investigated in this study and the properties of the corresponding datasets. 
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" See Sect. 12.1 I for the selection criteria and Sect.|4]for references concerning the source properties. 

* PWN = pulsar wind nebula, yuQ = microquasar, NS = binary containing a neutron star, BHC = containing a black hole candidate, ? = 
classification uncertain 

Calculated at 30 keV and taking into account the SPl response 

Assumed variability timescale in the analysis at £ < 144 keV (first number) and at £ > 144 keV (second number); "-" = the source was 
assumed constant either because it is known to be steady or because the variability timescale and/or amplitude is too small for SPI to detect 
significant variability. 

' Control position for the assessment of systematic errors in the analysis of crowded fields (see Sects. [33] and l4.2ll . 



put catalog in order to have a complete model. But as the en- 
ergy increases, many sources in the catalog drop below our flux 
sensitivity and need no longer be included. In order to optimise 
the sensitivity of our analysis, we therefore construct three sep- 
arate catalogs for three energy ranges: (a) 25-144 keV, (b) 144- 
346 keV, and (c) 346-8000 keV. For range (a), we use the en- 
tire Bouchet catalog. For range (b) we use only those sources 
from the Bouchet catalog which are detected at 100-200 keV 
(i.e. a flux is stated in BOS, not an upper limit) and which have 
at least 3 mCrab at 50-100 keV. Finally, for range (c) we include 
only those objects from the Bouchet catalog which have a flux > 
6 mCrab at 100-200 keV or which are marked as variable and de- 
tected at 50-100 keV. Figure[T]shows the positions of the sources 
in these three input catalogs in galactic coordinates. 

The fraction of the sources from the input catalogs which 
are actually included in the model constructed by spimodfit for 
an individual dataset, depends on how crowded the field is and 
how widely spread the pointings are within the search radius. In 
our datasets (see Table [U the number of sources included in the 
model ranges from 1 (in energy range (c) for the least crowded 
fields) to 83 (in energy range (a) of the most extended among 
our datasets, the one containing GRS 1758-258). 

3.1.3. Time variability 

The source variability is an important parameter in the flux ex- 
traction process. Bright and strongly variable sources can in- 
troduce some spurious features in the analysis if their model 
variability timescale is chosen too long and their true variabil- 



ity timescale and amplitude is large enough for SPI to detect a 
flux change. 

We treat those sources as time-variable which are marked as 
such in the B08 catalog. These are the following 12 sources (in 
brackets we give the variability timescale chosen by us based 
on the ISGRI light curves (Courvoisier et al. .2008) and the 
sensitivity of SPI): lA 0535+262 (180 ks), Vela X-1 (90 ks), 
GX 301-2 (90 ks), 4U 1700-377 (90 ks), IGR J17464-3213 
(90 ks), Sco X-1 (90 ks). Swift J1753.5-0127 (2.6 Ms), Aql X-1 
(180 ks), GRS 1915+105 (90 ks), Cyg X-1 (90 ks, 180 ks for 
E > 144 keV), EXO 0331 (2.6 Ms), Cyg X-3 (90 ks). Cyg X-1 
is the only source which is significantly variable over the entire 
energy range. All other variable sources are only treated as sig- 
nifcantly variable (given the SPI sensitivity) up to 144 keV. 

The minimum timescales used here were determined from 
a study of model fit convergence. We found that our model 
fits do not converge well if we assume timescales much below 
90 ks resulting in increased statistical eiTors on all fit parame- 
ters. Particularly in the case of the bright and strongly variable 
source 4U 1700-377, B08 used a variability timescale of »3 ks 
(the minimum possible due to the way the data is stored). For 
such sources, dedicated studies are necessary to further increase 
the accuracy of spectra for short integration times. However, it is 
our aim here to provide long-term average spectra with a coher- 
ent treatment of all sources. 

For the variable sources, the analysis yields a separate spec- 
trum for each variability time bin, or in other words a separate 
light curve for each energy bin. We derive the average spec- 
trum by applying a weighted least-squares procedure described 
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Fig. 1. (a) (top) The positions of the sources from the input catalogs used in this study for the analysis of the SPI data in three 
separate energy ranges: catalog 1 (25-144 keV, crosses, 173 objects, i.e. the entire BOS catalog), catalog 2 (144-346 keV, diamonds, 
129 objects), and catalog 3 (346-1000 keV, squares, 52 objects) on an all-sky map in galactic coordinates. The 20 objects investigated 
in this work are marked as thick squares. The grid mesh size 15°xl5°. Also shown are the extensions of the SPI analysis field of view 
and the ISGRI partially coded field of view, (b) (bottom) The enlarged galactic centre region of (a). The catalogs for higher energies 
are derived by successively removing weak sources from the initial catalog, i.e. a lower energy catalog contains all higher-energy 
ones. See Sect. 13.1.2] 



in Alvarez-Gaume et al. (I2004I I. page 14, to the light curves from 
each energy bin. 



factor equal to the number of energy bins, i.e. from several days 
to a few hours for each source data set. 



3.1.4. Postprocessing and spectral model fitting 

The analysis is sped up by processing each energy bin on a sep- 
arate processor reducing the required computing time by about a 



In a final step, the data from each processor are collected, 
where necessary the time-averaging of the spectra is performed, 
and the spectra are converted to the PHA format suitable for the 
spectral analysis with XSpec 12.3.1 (Arnaud et al. |2007ji which 
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is then used for the fitting of spectral models taking into account 
the SPI response and energy-redistribution. 

3.2. ISGRI data analysis 

The analysis of the INTEGRAL IBIS/ISGRI data is based on a 
cross-correlation procedure between the recorded image on the 
detector plane and a decoding array derived from the mask pat- 
tern (see Goldwurm et al. 2003 ). Standard spectral extraction has 
been applied as provided by the OSA 7.0 software package, us- 
ing the same energy binning and science windows as for the SPI 
data. 

3.3. Statistical and systematic errors 

The statistical errors of the photon count rates in each energy 
bin are determined both for the ISGRI and the SPI data by the fit 
software taking into account proper Poisson count statistics both 
for background and source signal. 

In addition there is a systematic eiTor on the flux calcula- 
tion stemming from the uncertainty of the calibration due to the 
limited amount of calibration data and the propagation of other 
uncertainties of the calibration process. We derive the relative 
error of the calibration of the spectra by assuming a shape for 
the Crab Nebula spectrum and independence of energy, i.e. that 
the uncertainty caused by the eiTor on the calibration is the same 
percentage of the measured count rate for each energy bin. This 
constant percentage is determined by varying it until the reduced 

from the fit of the assumed spectral shape to the Crab spec- 
trum is close to 1 .0 for the complete energy range. 

For both the ISGRI and the SPI data, we follow the 
INTEGRAL Cross-calibration status document by Jourdain et al. 
( I2008I I and assume the Crab Nebula spectrum between 25 keV 
and 1 MeV to be a broken powerlaw with break energy fixed at 
100 keV (which is close to the break energy we actually mea- 
sure, see Sect. 14. lb . With this assumption we derive a calibration 
error of our ISGRI and SPI spectra of 0.5 % in each energy bin. 
Unless otherwise noted, this error has been added in all spectral 
plots shown in this paper For the fainter sources, it is, however, 
negligible. 

The Crab is a bright source in a relatively quiet field without 
many other neighbouring hard X-ray sources. In order to assess 
the possible systematic eiTors arising in the analysis of crowded 
fields in the galactic bulge where some of our 20 sources of in- 
terest are located, we have added a control position to our input 
catalog at RA = 259.5°, Dec = -40.35° (1 = 347.7°, b = -2.1°). 
This position was randomly chosen in an empty region (with re- 
spect to the BOS catalog) of the crowded field near GX 354-0. 
The analysis for this artificial reference source is discussed in 
Sect ion 

4. Results 

In this section we present our spectral analysis results and dis- 
cuss the agreement of our SPI and ISGRI spectra with each other 
and with previously published measurements for each object in 
detail. The numerical results are later summarised in tables |2l [3] 
and|5]in Sect. |5] where we also discuss the agreement with the 
fluxes given by BOS. 

Studying data below 25 keV would exceed the scope of this 
paper. Also, due to the small field-of-view of the INTEGRAL 
soft X-ray instrument Jem-X, a soft-X-ray dataset truly concur- 
rent with our SPI and ISGRI dataset would have meant shrinking 



the SPI and ISGRI exposure by more than a factor two. When 
we compare to previously published measurements, this lack of 
soft X-ray data leads often to poorly constrained low energy fea- 
tures such as cutoffs near our lower energy limit (25 keV) or soft 
photon temperatures in thermal Comptonisation models. Where 
necessary, we try to overcome this by taking the appropriate val- 
ues from the literature. 

In the following, unless otherwise noted, the normalisation 
parameters "Norm" or "Normi" in the spectral fit results are 
given in units of photons keV 'cm^^s ^ They correspond to the 
flux at 1 keV. The letter F designates a powerlaw spectral in- 
dex. If we need to introduce a second power-law component for 
high-energy emission, we use the XSpec model PEGPWRLW 
with peg energy 200 keV. The normalization of it ("Norm2") is 
given as the differential flux at 200 keV. 

In the figures in this section we show the spectra unfolded 
by XSpec with statistical eiTors (slightly increased by the cali- 
bration errors discussed in Sect. l3.3b . The solid line always indi- 
cates the fitted model (as a step function according to the energy 
binning). If the model has several components, these are also 
separately shown with dashed lines. 

4.1. Crab Nebula 

The Crab Nebula is the standard candle of high-energy astron- 
omy. Because of the steady, strong emission from this pulsar 
wind nebula, it has always been among the first objects to be 
observed by new instruments. For a recent compilation of X-ray 
data on the Crab see, e.g., Kirsch et al. (l2005b . Still, its high en- 
ergy spectrum is not too well known. Many instruments use an 
assumption about the Crab's spectral shape as the basis of their 
calibration and can hence not say much about the fine details of 
the true spectrum. This is also the case for INTEGRALflSGRl. 
On the other hand, INTEGRAL/SPl with its high energy resolu- 
tion could measure a Crab spectrum which is less dependent on 
such assumptions about the spectral shape. 

Figure |2^ shows the result of our analysis of the SPI Crab 
spectrum. A good fit is a broken powerlaw (XSpec model "bkn- 
pow") with the following parameters: 
Fi = 2.11 ± 0.01 

£break = (SI ± ll)keV 

F2 = 2.20 + 0.01 
Norm = 11.6 + 0.5 

with a reduced of O.S for 13 degrees of freedom (see also 
Sect. [331). 

The corresponding ISGRI dataset gives the following Crab 
spectrum (Fig.Eb): F, = 2.1 19 ± 0.006, £break = (96+ 4) keV, F2 
= 2.36 + 0.02 , Norm = 10.7 + 0.3, with a reduced of 1.3 for 
11 d.o.f. (see Sect.O. 

This needs to be compared to our best fit for the SPI data. 
The shape is in very good agreement below 100 keV. Table |4] 
(see Sect. |5]) gives our measurements of the Crab flux in sev- 
eral energy bands. The model fluxes between 25 and 100 keV 
( 0.22S7±0.0004 cm^^g-i for SPI and 0.2046+0.0004 cm^^s-' 
for ISGRI ) diff^er by a factor 1.1 1+0.01 (SPI/ISGRI) which has 
to be regarded as a measurement of the systematic difference in 
absolute flux normalisation. 

Above the fixed break, the ISGRI spectrum is significantly 
softer than the SPI spectrum. However, as the ISGRI response 
was made (by the authors of the OSA 7 ISGRI response) to 
give the Crab spectrum a powerlaw above 100 keV, this devi- 
ation tells us that there are probably significant systematic errors 
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Fig. 2. The unfolded spectra of the Crab Nebula as derived from 
dataset 1 (see table [TJ with statistical errors: (a) (top) SPI spec- 
trum (b) (bottom) ISGRI spectrum; fit function is in both cases a 
broken powerlaw. Parameters of both fits are given in Sect. 14.11 



Fig. 3. The unfolded spectra of the Vela Pulsar and Pulsar Wind 
Nebula derived from dataset 2 (see table[T]l with statistical errors: 
(a) (top) SPI spectrum, (b) (bottom) ISGRI spectrum. Fit func- 
tion is in both cases a simple powerlaw (parameters are given in 
Sect.O- 



in the ISGRI calibration at the upper end of its energy range in 
addition to the calibration errors derived by us in Sect. 13.31 

Our SPI spectrum is in very good agreement with that de- 
rived from SPI data up to 1 MeV by Jourdain et al. (2008) whose 
analysis was not based on spimodfit. Furthermore, the smooth fit 
of the narrow energy bin around the 511 ke V line, which is also a 
strong background line of the instrument, demonstrates that the 
background determination works well. 

4.2. Vela Pulsar 

The other pulsar wind nebula/pulsar in our set of sources, 
the Vela PWN, is generally regarded as similar to the Crab. 
However, in hard X-rays the object is about two orders of mag- 
nitude fainter even though it is nearly an order of magnitude 
closer It has been studied in detail by all high-resolution soft 
X-ray observatories (see e.g. Mangano et al. (2006) and refer- 
ences therein). 

Figure [3^ shows the spectrum from our analysis of the SPI 
Vela Pulsar data. Because of the weakness of the source, we 
choose the wider energy binning in the SPI analysis. A simple 



powerlaw with F = 2.13 +0.10, Norm - 0.12+0.05 describes the 
spectrum well. The reduced of the fit is 0.55 for 9 d.o.f. 

The SPI spectrum is in good agreement with our correspond- 
ing ISGRI spectrum (Fig. |3j)) which has the parameters F = 
2.04+0.04 , Norm = (0.057+0.008), with a reduced of 1.4 
for 13 d.o.f. 

The hard X-ray spectrum of this source was studied up to 
200 keV by Mangano et al. (2006) using Beppo-SAX data. They 
find a spectral index of 2.00+0.05 at 15-200 keV which is in very 
good agreement with our measurements. 

At energies above 200 keV, our analysis does not show sig- 
nificant evidence for emission. On the other hand, our measure- 
ment is well consistent with a continuation of the spectrum in a 
way similar to what we observe for the Crab. Given the nature 
of the source and its morphological similarity to the Crab, this is 
to be expected. OSSE detected the Vela PWN only marginally 
between 200 and 760 keV with a spectral index of 1.8+0.3 
(Strickman et al. |1996bl ) compatible with a continuation of our 
spectrum without a break. 
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Fig. 4. The unfolded spectra of NGC 4151 derived from dataset 
3 (which has very low exposure, see table [U with statistical er- 
rors: (a) (top) SPI spectrum, (b) (bottom) ISGRI spectrum. The 
fit function is in both cases a powerlaw (parameters are given in 
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Fig. 5. The unfolded spectra of NGC 4945 derived from dataset 
4 (which has low exposure, see table[T]l with statistical errors: (a) 
(top) SPI spectrum, (b) (bottom) ISGRI spectrum. The fit func- 
tion is in both cases a powerlaw with exponential cutoff and fixed 
photo-electric absorption (parameters are given in Sect. 14.4b . 



4.3. NGC 4151 

The prototypical Seyfert galaxy NGC 4151 is one of the best 
studied active galaxies. For a comprehensive review on this ob- 
ject see UMch (2000). Unfortunately, this is our smallest dataset 
with less than 100 ks observation time. We show the result- 
ing spectra here only for completeness. A detailed study of this 
source using INTEGRAL data up to 300 keV was published by 
Beckmann et al. (2005) (they use only ISGRI data and get a 
larger dataset by using data taken during SPI annealings). 

Figure |4^ shows the SPI spectrum. The best fit is a simple 
powerlaw with F = 2.10+0.17, Norm - 0.38+0.25, and a reduced 

of 1.4 for 9 d.o.f. The excess in the annihilation line bin is 
marginal (2.4 cr significance). 

The ISGRI spectrum (Fig. |4j5) gives a compatible fit result 
with r = 2.03+0.06, Norm = 0.24+0.05, reduced ;r^=1.4 for 13 
d.o.f. 

Fitting a powerlaw with exponential cutoff to the ISGRI data 
as in Beckmann et al. (2005) improves the fit slightly (x^-1.3) 
and results in values well compatible with previously published 



measurements: F — 1.57+0.24, -Ecutoff = 121+55 keV, Norm = 
0.07+0.04 . 

4.4. NGC 4945 

NGC 4945 is the brightest Seyfert 2 and the second brightest 
hard X-ray source (after NGC 4151) of all radio-quiet AGN 
(Done et al. 1996 ). Its X-ray spectrum has been extensively stud- 
ied by various authors (see e.g. Itoh et al. (2008 ) and references 
therein). As a Seyfert 2, the source is strongly absorbed at lower 
energies. Itoh et al. (120081 ) using Suzaku measure a hydrogen 
column density A^h of 5.4x 10^"* cm~^ and a hard X-ray spectrum 
with a powerlaw index of 1.5+g j and an exponential cutoff at 
an energy of 150+5™' ^™ agreement with previous measure- 
ments). 

Our dataset for NGC4945, like the one for NGC 4151, has 
relatively low exposure. The hydrogen column density for this 
source is by far the highest of all the sources in our sample. We 
find that it is still low enough to be negligible for our SPI spec- 
tral modelling but the ISGRI sensitivity below 50 keV is high 
enough to make the spectral fits sensitive to the presence of the 
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Fig. 6. The unfolded spectra of Cen A derived from dataset 5 
(see table [U with statistical errors: (a) (top) SPI spectrum, (b) 
(bottom) ISGRI spectrum. The fit function is in both cases a 
powerlaw (parameters are given in Sect. 14.5b . 
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Fig. 7. The unfolded spectra of XTE J 1550-564 derived from 
dataset 6 (see table [T]) with statistical errors: (a) (top) SPI spec- 
trum, (b) (bottom) ISGRI spectrum (increased errors were re- 
quired, see text). The fit function is in both cases a powerlaw 
with exponential cutoff (parameters are given in Sect. l4.6l l. 



absorption. In order to compare better with previous measure- 
ments, we therefore include the photoelectric absorption in the 
spectral model (Xspec model WABS) and fit an absorbed pow- 
erlaw with exponential cutoff to the SPI and the ISGRI spectrum 
fixing the hydrogen column and the cutoff energy to the values 
found by Itoh et al. (|2008|) . 

Figure |5h shows the resulting SPI spectrum. The best fit pa- 
rameters are F = 1.66+0.17, Norm = 0.061+0.039 with a re- 
duced of 0.4 for 9 d.o.f. 

The ISGRI spectrum is shown in Fig. [S]?. Keeping A^h fixed 
as for SPI, the best fit parameters are F - 1.56+0.05, Norm - 
0.033+0.006 with a reduced of 0.6 for 13 d.o.f. 

SPI and ISGRI spectrum are both consistent with each other 
and with the spectra measured by Itoh et al. (12008l l. 

4.5. Cen A 

The nearby radio galaxy Cen A has been observed by all X- 
ray observatories (see Rothschild et al. (2006) and references 
therein). It has been classified as a Seyfert 2 (see Dermer & 
Gehrels ri995l l. The consistently observed hard X-ray spectrum 
is, in spite of flux variability, a simple powerlaw of index w 1.8 



(Rothschild et al. |2006l which extends to the MeV range where 
the spectrum steepens (Steinle et al. ll998l l. 

Even though the available SPI exposure for Cen A is mod- 
est, we obtain a clear detection also beyond 200 keV. The SPI 
spectrum (Fig.|6h) gives the best fit parameters F - 1.89+0.05, 
Norm = 0.20+0.04, and a reduced of 1.2 for 9 d.o.f.. 

Similarly, the ISGRI spectrum (Fig.jSJ?) gives a best fit with 
F = 1.80+0.02, Norm = 0.1 17+0.009 and a reduced of 0.9 for 
13 d.o.f. 

SPI and ISGRI spectral shape are in excellent agreement 
with each other and previous measurements. 

4.6. XTE J 1550-564 

The Low-Mass X-ray Binary XTE J 1550-564 has been classi- 
fied as a microquasar and a black hole candidate. Since its dis- 
covery in 1998 (Smith 1998) it has been subject of several de- 
tailed studies and was found to be similar to Cyg X-1, however, 
with lower luminosity (see e.g. Wu et al. 120071 and references 
therein). Sturner & Shrader ( I2005I I studied the source using 
RXTE and INTEGRAL data and found that it stayed in its low- 
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hard spectral state even when going through modest outbursts. 
According to the ISGRI lightcurve of the source (Courvoisier 
et al. 12008 1). XTE J1550-564 had only two short major outbursts 
with INTEGRAL coverage and was otherwise in a low state. The 
time-averaged spectrum we are determining here should there- 
fore be close to the low-hard type, i.e. consistent with thermal 
Comptonisation which is essentially a powerlaw with exponen- 
tial cutoff in our energy range. 

We find that a powerlaw with exponential cutoff fits the SPI 
data reasonably. From the SPI spectrum (Fig. |7^) we get the 

parameters F = 1.36+0.09, ^cutoff = (206±50) keV, Norm = 
0.031+0.010 with a reduced of j 5 ^j. g ^ f 

For the ISGRI spectrum we don't achieve a reduced be- 
low 2.0 with any common model (we tried cutoff powerlaw, ther- 
mal Comptonisation, and black body + powerlaw) although the 
residuals don't show any clear trend in the deviations. This can 
be explained by additional systematic errors in the image recon- 
struction caused by the complex source region. Allowing for this 
by doubling the assumed calibration eiTors from 0.5 % to 1.0 % 
(see also Sturner & Shrader (I2005I I). we get a reduced of 1.3 
for 12 d.o.f. The best fit parameters are marginally consistent 
with the ones obtained from the SPI data: F - 1.00+0.05, E^^oS 
= (108+9) keV, Norm = 0.009+0.002. This is due to higher pho- 
ton fluxes measured by ISGRI below 200 keV. Above 200 keV, 
the fluxes are consistent (cf. Table|5]). 

For comparison with earlier publications, we also fit the 
thermal Comptonisation model (COMPTT) by Titarchuk ( |1994| l 
which is the best fit to the ISGRI and SPI data in Sturner 
& Shrader (2005). Averaging over their datasets, they obtain 
a plasma temperature kT - (48.4+2.7) keV, a plasma optical 
depth Tp - 1.48+0.07, and a soft photon temperature kTo - 
(0.5+0.14) keV. Since we lack the data between 3 keV and 
25 keV, kTo is not well constrained. We therefore fix it to the 
value 0.5 keV by Sturner & Shrader. From the SPI data we 
then obtain kT = (56.9+7.7) keV and Tp = 1.32+0.18 (reduced 
X^-13 for 9 d.o.f.) in good agreement with Sturner & Shrader 

4.7. 4U 1630-47 

The X-ray transient 4U 1630-47 has been classified as a LMXB, 
a very likely black hole candidate, and has also been found to 
show indications of the presence of jets but is not yet regarded as 
a microquasar (see Kubota et al. (2007) and references therein). 
During the observation time covered by our dataset, the source 
showed frequent but short outbursts. The time-averaged spec- 
trum should therefore be a superposition of the various states the 
source can be in, with the low state dominating. Tomsick et al. 
(l2005T l analyse RXTE and INTEGRAL data from 4U 1630-47. 
They distinguish five different spectral states which they all fit 
with a two-component model made from a black body with inter- 
stellar absorption and a powerlaw. They find that the parameters 
of this model vary strongly between the states. The lowest flux 
state is found to be the one where the emission is dominated by 
black body emission {kT a; 1.4 keV) below 15 keV, and above 
that by a steep powerlaw with index F ^3. A simple powerlaw is 
also the best fit to our spectra of this on average relatively weak 
source. 

The SPI spectrum for 4U 1630-47 (Fig. [D is best fit by a 
powerlaw with F = 3.3+0.2 and Norm - 33+21. The reduced;^^^ 
is 0.6 for 9 d.o.f. 

The ISGRI spectrum may suffer from source contamination 
from the neighbouring source IGR J16358-4726, especially at 
higher energies (see also Tomsick et al. I20051 l. We try to ac- 
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Fig. 8. The unfolded spectrum of 4U 1630-47 derived from the 
SPI dataset 7 (see table[Tl) with statistical errors. The correspond- 
ing ISGRI data is flawed because of source contamination by 
IGR J16358-472 (see text). The fit function is a powerlaw (pa- 
rameters are given in Sect. l4.7l) . 

count for this approximately by increasing the systematic errors 
to 5 %. Fitting a powerlaw then results in F = 2.33+0.24, Norm 
= 0.04+0.04 with a reduced = 1.0 for 13 d.o.f. which is much 
harder than the SPI spectrum. 

Furthermore, the observed fluxes (see Table|5]l from SPI and 
ISGRI differ strongly at low energies while the SPI measure- 
ments agree with those from B08 who also note that the source 
may be confused with IGR J16358-472. 

We conclude that the neighbouring source IGR J16358-472 
may be hampering the measurement for ISGRI or SPI or both. 
Since this distortion is strongly dependent on the instrument 
characteristics, the effects are different for SPI and ISGRI and 
lead very probably to the observed difference in the spectra. The 
smaller differences between our spectra and those from Bouchet 
are to be expected since here the same instrument is used. The 
fact that the SPI spectrum agrees well with RXTE measurements 
of 4U 1630-47 in a low state by Tomsick et al. (2005) argues for 
the better quality of the SPI spectrum. Also, our source model 
for the SPI data explicitely includes IGR J16358-472. The SPI 
analysis should therefore in this case be more immune to con- 
tamination than the ISGRI analysis. 

4.8. Swift J1 656.3-3302 

The recently discovered source Swift J1656. 3-3302 (Okajima et 
al. l2006l l has been identified with a high redshift (z = 2.4) blazar 
(Masetti et al 12008) 1. The late discovery of this source is due to 
its hard spectrum, i.e. its weakness at soft X-rays, and to its loca- 
tion near the galactic bulge which makes for a very complex field 
of view with many variable and bright neighbouring sources. 
Probably because many bright sources are very near the edge 
of the partially coded field of view when Swift J1656. 3-3302 is 
near its centre, our standard modelling technique turned out to 
be inadequate for extracting a SPI spectrum of Swift J1656.3- 
3302 from a dataset centred at the source. The quality of the spi- 
modfit fits was never acceptable in more than a few energy bins. 
But when the source was observed by SPI with high exposure 
and off-axis by at least 5° in the dataset selected for GRS 1758- 
258 (dataset 10, see table[T]), the extraction of a useful spectrum 
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Fig. 9. The unfolded SPI spectrum of Swift J1656.3-3302 de- 
rived from dataset 10 (see table[T]i with statistical errors. The fit 
function is a powerlaw (parameters are given in Sect. l4.8l l. 



with our standard input catalog and variability information was 
possible in all energy bins except for the lowest two. The cor- 
responding ISGRI spectrum of Swift J1656. 3-3302 from dataset 
10, however, was not useful since the object was too far outside 
the ISGRI fully coded field of view. We therefore present in Fig. 
|9]only the SPI spectrum for this source. It is well described by a 
simple powerlaw with T - 1.49+0.18 and Norm = 0.004+0.005. 
The reduced is 1.0 for 7 d.o.f. 

An ISGRI spectrum for Swift J1656. 3-3302 was measured 
by Masetti et al. (2008 ). They find a spectral index of 1.64+0.16 
in agreement with our measurement. 
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4.9. OAO 1657-415 

The high mass X-ray binary OAO 1657-415 is a system consist- 
ing of an eclipsing 38 s pulsar and a highly reddened B super- 
giant with orbital period 10.4 d (see e.g. Audley et al. '2006' and 
references therein). The object has been studied since its discov- 
ery in the 1970s by a fair number of authors with several X-ray 
observatories up to about 160 keV, recently also with INTEGRAL 
(Filippova et al. 120051 Barnstedt et al. l2008l) . In the range from 
3 to 160 keV they find consistently an absorbed, exceptionally 
hard spectrum with index F ^ 1.0 which then softens with an 
exponential cutofi" at energies around 25 keV. 

We fit our time-averaged (and therefore also phase-averaged) 
SPI spectrum initially with the simplest of these models, the 
powerlaw with exponential cutoff. We obtain F = 1 .17+0.23 and 
^cutoff = (23.9+3.2) keV with a reduced of 1.4 for 8 d.o.f., in 
good agreement with previous measurements. 

Inspecting the residuals, we find marginal evidence for a 
hard tail diverging from the model above 150 keV. The feature 
is not excluded by previous studies since they did not cover this 
energy range. The significance of the hard tail is about 3 cr which 
does not warrant a much more detailed study. In order to quan- 
tify the possible flux of this feature, we add a powerlaw to the 
model keeping E cutoff fixed. This results in a lower reduced of 
0.7 and best fit pai-ameters Fi = 1.23+0.07, Norm i =0.16+0.04, 
F2 = 1.27+0.58 and a flux of the second component at 200 keV 
of (2.1 + 1.0) X 10-^ cm-^s-ikeV-i. Fig. [Toll shows the fit result 
with the added powerlaw. 



Fig. 10. The unfolded specti-a of OAO 1657-415 derived from 
dataset 7 (see table [T]) with statistical errors: (a) (top) SPI spec- 
trum fit with powerlaw with exponential cutoff and additional 
powerlaw component (b) (bottom) ISGRI spectrum fit with pow- 
erlaw with exponential cutoff (with increased errors, see text). 
The fit parameters are given in Sect. 14.91 

For the ISGRI data, we increase the statistical errors by 5 % 
of the flux in order to get take into account the larger systematic 
uncertainties of the analysis (this is the same dataset as the one 
for 4U 1630-47 where this increase was also necessary). 

If we fit the same model with added powerlaw as for the SPI 
data, we get reduced of 0.4 for 10 d.o.f. compatible with the 
SPI fit. The parameters of the second powerlaw are essentially 
unconstrained. The reality of the potential hard tail in the SPI 
data is not ruled out. 

If we fit with the powerlaw with exponential cutoff without 
the second powerlaw, we get a reduced of 0.4 for 12 d.o.f. 
and best fit parameters Fi = 1.4+0.3, /icutoff = 25+4 keV, Norm 
- 0.23+0.21 which is well compatible with the SPI data. This 
spectrum in shown in Fig.llOb. 

4.10. GX 339-4 

The microquasar GX 339-4 in one of the classical black hole X- 
ray binaries and has been studied by all X-ray and soft gamma- 
ray observatories, also with Swift (Tomsick et al. 1200 8I 1 and 
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Fig. 11. The unfolded spectra of GX 339-4 derived from dataset 
8 (see table [U with statistical errors: (a) (middle) SPI spectrum 
fit with two-component model: powerlaw with exponential cut- 
off and second powerlaw. (b) (bottom) ISGRI spectrum fit with 
the same model as in (a). The fit parameters are given in Sect. 

INTEGRAL (see Joinet et al. l2007l who also give a review of 
the literature). Like other black hole binaries, the source has 
two main states, the low-hard state during which the high-energy 
spectrum is essentially a powerlaw with index F = 1.4-2.1 with 
an exponential cutoff at hundred to a few hundred keV, and the 
high-soft state during which T goes to values > 2.4 but without 
cutoff (see again Joinet et al. '2007). Studies of the soft-gamma 
ray emission with OSSE have also found emission up to 400 keV 
(Grabelsky et al. I1995I I. Joinet et al. (I2007I I found evidence for 
emission above 200 keV in the low-hard state at the 5 cr level. 

Our time-averaged spectrum should be a superposition of all 
states of the object. As a first step, we fit a simple powerlaw 
to the data and obtain index T = 2.08+0.04, Norm = 0.35+0.05 
with a reduced of 1.2 for 15 d.o.f. The residuals are systemat- 
ically negative from 100-300 keV but seem to be systematically 
positive above 300 keV indicating a hard tail. The significance 
of the emission above 200 keV is 6.4cr. 

The same fit applied to the ISGRI data gives a good reduced 
of 0.6 for 13 d.o.f and best fit parameters F - 2.07+0.04 and 
Norm = 0.10+0.01 confirming the shape of the SPI result up 



to 200 keV. Beyond 200 keV, the ISGRI result is not sensitive 
enough. As in the SPI spectrum, there is no evidence for a cutoff. 
A hard tail cannot be ruled out by the ISGRI data. 

In order to improve our fit to the SPI data and to better be able 
to compare with previous measurements, we fit a more complex 
model consisting of two components: (a) a powerlaw with ex- 
ponential cutoff (in our energy range this is equivalent to a ther- 
mal Comptonisation model, one of the main candidate emission 
mechanisms for microquasars, see Sect. 14.6b and (b) an addi- 
tional powerlaw in order to describe a potential hard tail which 
may obscure the observation of a cutoff in our spectrum. 

Since we lack the soft X-ray data to better constrain the 
cutoff energy, we fix it to 100 keV consistent with previous 
measurements (if iScutoff is permitted to vary freely, we get a 
best fit value of 88+56 keV). Figure [TTh shows the result of 
this fit. The best fit parameters are Fi = 1.6+0.1, Normi = 
0.09+0.03, F2 = 1.0+0.6, Norma (Flux at 200 keV) = (3.8 + 
2.4) X 10"'' cm"-s"'keV"' with an improved reduced of 0.8 
for 13 d.o.f. 

The index Fi of the dominant powerlaw at lower energies is 
well consistent with that observed by Tomsick et al. (2008) up to 
100 keV. The additional powerlaw becomes dominant at about 
200 keV. This confirms the feature akeady noted by Joinet et al. 
(I2007I I who speculated that it could be produced in a jet. 

The more complex model is also fully compatible with the 
ISGRI data. Fitting the two-component model with fixed iicutoff 
- 100 keV and Y2-I.Q to the ISGRI data, we obtain a reduced 

of 0.7 for 12 d.o.f. and best fit parameters Fi - 1.7+0.1 (Fig. 

4.11. 4U 1700-377 

The eclipsing, 3.4 day X-ray binary 4U 1700-377 consists of a 
hot, bright supergiant and very probably a neutron star because 
its overall X-ray spectrum is similar to that of an accreting pul- 
sar The system has been studied in detail by many authors. For 
a review see e.g. van der Meer et al. ( 120051) . In soft X-rays, 
the source displays a multitude of line features. The high en- 
ergy spectrum has been modelled successfully by a powerlaw 
with high energy cutoflF with additional photoelectric absorption 
for the points below ^ 25 keV. Authors have also employed 
more complex models like thermal bremsstrahlung or thermal 
Comptonisation. In all cases, even though the source is strongly 
variable, the spectral shape above about 25 keV was found to be 
independent of the flux state. 

Our SPI spectrum is indeed well described by a powerlaw 
with high energy cutoff (Xspec model HIGHECUT). Figure [T2k 
shows the SPI spectrum. The best fit parameters are F = 2.3+0.1, 
Norm = 5.0+2.4, E,. = (31+2) keV, £/ = (65+11) keV. The re- 
duced is 0.9 for 7 d.o.f. 

This agrees well with the literature: e.g. Maisack et al. ( 19^ 
find F = 2.55!° '^, E^ = (20. 1+'' ;|) keV, £/ = (69.8+1-^^) keV. But 
also a simpler power law with exponential cutoff at 48+4 keV 
and F = 2.0+0. 1 describes the data reasonably. 

As for the other two sources in dataset 7, the statistical errors 
for the ISGRI spectrum of 4U 1700-377 (Fig. [Eb) had to be 
increased by 5 % of the flux. This is certainly at least partially 
caused by the strong variability of 4U 1700-377 (see also Sect. 
13.1.31 ). A fit with a power law and exponential cutoff gives the 
best fit parameter values F - 2.3+0.2, Norm = 7+3, -Ecutoff = 
(74+14) keV, in agreement with SPI. The reduced is 0.9 for 
12 d.o.f. 



Petry, Beckmann, Halloin, Strong: Soft gamma-ray sources detected by INTEGRAL 



13 



4U 1700-377 (SPI) 



IGRJ1 7091 -3624 (SPI) 




100 200 
Energy (keV) 



4U 1700-377 (ISGRI) 



100 200 
Energy (keV) 




100 200 
Energy (keV) 



> IGR Jl 7091 -3624 (ISGRI) 




100 200 
Energy (keV) 



Fig. 12. The unfolded spectra of 4U 1700-377 derived from 
dataset 7 (see table [U with statistical errors: (a) (top) SPI spec- 
trum, (b) (bottom) ISGRI spectrum. Spectral model is in both 
cases a powerlaw with high energy cutoff. The fit parameters are 
given in Sect. 14.111 



Fig. 13. The unfolded spectra of IGR J 1 709 1 -3624 derived from 
dataset 9 (see table [T]) with statistical errors: (a) (top) SPI spec- 
trum with powerlaw fit, (b) (middle) ISGRI spectrum with pow- 
erlaw fit. The fit parameters are given in Sect. 14.121 



4.12. IGR J1 7091 -3624 

The object IGR J17091-3624 was discovered in 2003 by 
INTEGRAL (Kuulkers et al. 2003) and subsequently studied by 
a small number of authors based on archival data (Revnivtsev 
et al l2003bT) and I NTEG RAL and KXTE data (Capitanio et al. 
I2006I I. Chaty et al. (I2008I I recendy found that the most probable 
candidate for a companion star is neither a giant nor a super- 
giant. They conclude that the system is most probably a LMXB. 
The source is variable but weak (< 10 mCrab above 25 keV). 
The spectrum above a few ten keV is in good approximation a 
hard powerlaw with index varying between 1.6 and 2.2, but has 
also been successfully modelled with thermal Comptonisation 
(Capitanio et al.l2006i 

When fitting our time-averaged SPI spectrum of 
IGR J17091-3624 (Fig. [T3li), we find that a simple power- 
law describes it well. We find as best fit parameters F = 1.8±0.2, 
Norm = 0.02+0.01 with a reduced x"^ of 0.4 for 9 d.o.f. This 
is consistent with the previous measurements mentioned above 
and with the notion that the source was mostly in the low-hard 
state during our observation time. 



Fitting a powerlaw to the corresponding ISGRI spectrum 
(Fig. [131)) gives a reduced of 1.2 for 16 d.o.f.,r = 1.93+0.04, 
and Norm = 0.034+0.005 consistent with the SPI spectrum. 

In order to compare our spectra with the ISGRI spectra by 
Captanio et al. (2006), we also fit a thermal Comptonisation 
model (Xspec model COMPTT). We follow Captanio et al. 
and fix the soft photon temperature ^T^o = 0.1 keV. From our 
ISGRI spectrum we then obtain a plasma optical depth Tp - 
1.5+0.3, and a plasma temperature kT - (36+8) keV, Norm 
= 0.0016+0.0005, and a reduced of 0.5 for 12 d.o.f. (Fig. 
[T3b). This is consistent with the values found by Captanio et al. 
(I2056l) . 

4.13. GX 354-0 

The LMXB GX 354-0, also known as 4U 1728-34, consists of an 
accreting neutron star and probably a main sequence star and has 
been classified as a bursting atoll source which implies signifi- 
cant spectral variability. It has been studied by most X-ray obser- 
vatories including INTEGRAL (see Falanga et al. 120061 and ref- 
erences therein). Above 25 keV, the spectrum has been typically 



14 



GX 354-0 (SPI) 



Petry, Beckmann, Halloin, Strong: Soft gamma-ray sources detected by INTEGRAL 

„ 1E 1740.7-2942 (SPI) 



100 200 
Energy (keV) 



GX 354-0 (ISGRI) 




100 200 
Energy (keV) 



Fig. 14. The unfolded spectra of GX 354-0 derived from dataset 
9 (see table[T]l with statistical: (a) (top) SPI spectrum with pow- 
erlaw fit, (b) (bottom) ISGRI spectrum fit with powerlaw. The fit 
parameters are given in Sect. 14.131 The combined significance of 
the emission above 200 keV from the SPI and ISGRI data is 3 cr. 
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Fig. 15. The unfolded spectra of IE 1740.7-2942 derived from 
dataset 9 (see table [T]) with statistical errors: (a) (top) SPI spec- 
trum, (b) (bottom) ISGRI spectrum, in both cases fit by a power- 
law with exponential cutoff. The fit parameters are given in Sect. 

KM 



modelled with a thermal Comptonisation spectrum. However, 
the studies do not extend much beyond 100 keV. 

In our spectral analysis of the time-averaged spectrum, we 
initially fit a simple powerlaw. For the SPI data, this already 
gives a satisfactory fit (Fig. fT4h ) with a reduced of 0.7 for 
9 d.o.f. and best fit parameters F = 3.26+0.06 and Norm - 
26.4+5.9 . The same is true if we fit a powerlaw to the ISGRI 
data (Fig.fT4b): We obtain a reduced of 1.0 for 13 d.o.f., and 
parameters F - 3.53+0.02, Norm = 43.3+3.0, steeper than the 
SPI spectrum. 

Falanga et al.C2006) find that a thermal Comptonisation 
model (Xspec model COMPTT) describes the data well. 
Lacking soft X-ray information we take the values for the soft 
photon temperature kT^) and the plasma temperature kT from the 
recent measurement by Falanga et al. (for the low-hard state), 
^ro=1.18 keV and kT =35 keV, and fix them in our fit. Applied 
to the SPI data this results in a reduced of 1.9 for 9 d.o.f. and 
unevenly scattered residuals suggesting a hard tail. However, the 
sensitivity of our measurement does not permit a more detailed 
study of this feature. Also the ISGRI data are not compatible 
with the pure COMPTT spectrum proposed by Falanga et al. 



Here we get a reduced of 7.6 for 11 d.o.f. even if we let all 
model parameters vary freely. 

We conclude that our data are not well compatible with a 
pure thermal Comptonisation spectrum since it underpredicts the 
high-energy emission. The combined significance of the emis- 
sion above 200 keV from GX 354-0 is, however, only 3 cr. 
Possibly, the averaging over several spectral states mimics a 
power law. 

4.14. 1E 1740.7-2942 

The black hole candidate IE 1740.7-2942 is a bright persis- 
tent source very near the Galactic Centre and the first source 
to be classified as a microquasar (see Bosch-Ramon et al. 120061 
for a recent review). Its classification as a LMXB is controver- 
sial. Like other BHCs, IE 1740.7-2942 exhibits low-hard and 
high-soft states. Using ISGRI data, Del Santo et al. (2005) find 
that the average spectrum above 25 keV is approximately de- 
scribed by a powerlaw of index F = 2.3 . In individual peri- 
ods they find harder spectra (F » 1.5) with exponential cutoffs 
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around 100 keV which can alternatively be described with ther- 
mal Comptonisation models. 

Bouchet et al. ( 1991), using G/MA^Ar/SlGMA data, have re- 
ported a transient additional component (modelled by a Gaussian 
with FWHM of 240 keV) peaking at 480 keV which appeared for 
about a day in the spectrum and then disappeared. This observa- 
tion has not yet been confirmed and our study would not pick 
up on such features since we are averaging over much longer 
timescales (unless such events would occur frequently). 

We fit our SPl spectrum with a powerlaw and exponential 
cutoff (Fig.fTSh) and obtain a reduced of 1.4 for 10 d.o.f. with 
r = 2.17+0.14 in agreement with the average spectrum by Del 
Santo et al. (2005 ). The normaHsation is Norm = 0.9+0.4. The 
cutoff energy is a badly constrained 434+400 keV. 

The same model fit to the ISGRI data gives best fit pa- 
rameters of r = 1.44+0.04, ^cutoff = (113+8) keV, Norm = 
0.040+0.005, and a reduced of 0.7 for 12 d.o.f. 

The ISGRI spectrum is in agreement with spectra found by 
Del Santo et al. (2005 ) (see above) but the agreement with SPI 
is at best marginal. In order to reconcile the two results, we fix 
the exponential cutoff in the model for the SPl data at the value 
determined by ISGRI. We obtain a reduced of 1.9 for 9 d.o.f. 
and best fit parameters F = 1.80+0.06, Norm = 0.28+0.06. The 
residuals suggest an additional hard component. But our sensi- 
tivity is not good enough to warrant a more detailed study. The 
agreement of the SPl results with the fluxes reported by B08 is 
very good (see Table|5]l. 

4.15. IGRJ1 7464-321 3 

The transient IGR J 1746 4-32 13 was discovered by INTEGRAL 
in 2003 (Revnivtsev et al. l2003al l and soon identified as the coun- 
terpart of the HEAO 1 X-ray source H 1743-322. The source was 
subsequently studied by several authors at X-ray, optical and ra- 
dio wavelengths, in particular with INTEGRAL (see Joinet et al. 
12005 for a recent review). Evidence for jets was found and the 
source is now regarded to be a microquasar (Corbel et al. 2005). 
Timing and spectral properties of the object suggest that it is a 
LMXB containing a black hole of about 10 Mq. 

The RXTEIASM Ughtcurve of IGR J17464-3213 (e.g. 
Capitanio et al. 2006a) shows that the source can be in strong 
outburst for months at a time. In between outbursts it returns 
to a very low level of emission and remains there for months. 
Using data from RXTE and INTEGRAL up to 150 keV, Joinet 
et al. ( 120051) have studied the complex spectral variability of the 
object (typical of BHCs): the source exhibits softening during 
flux increases followed by phases of constant spectrum in spite 
of intensity changes by factors 2-3. Above a few ten keV, the 
spectrum is described by a powerlaw with index F varying be- 
tween 1.0 and 2.6, and a not very well constrained exponential 
cutoff around 100 keV. 

The ISGRI lightcurve (Courvoisier et al. 12008] ) for 
IGR J 17464-32 13 shows three strong outburst periods and oth- 
erwise a low level of emission such that our time-averaged spec- 
trum should be a superposition of the low-hard and the high-soft 
state. When examining our SPI spectrum, we find evidence for 
increased uncertainties in the lowest two energy bins possibly 
due to to insufficient source region modelling at the lowest en- 
ergies. To be on the safe side, we increase the statistical errors 
by 5 % of the flux for the entire spectrum. Fitting the spectrum 
with a powerlaw (Fig.[T6b). then leads to a reduced of 1.1 for 
9 d.o.f. and best fit parameters F = 2.82+0. 14, Norm = 2.1 + 1.1. 

Fitting a powerlaw to the ISGRI spectrum of (Fig.[T6t) gives 
a perfect reduced of 0.9 and best fit parameters F = 2.49+0.02 
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Fig. 16. The unfolded specti-a of IGR J17464-3213 (a.k.a. 
H 1743-322) derived from dataset 9 (see table[Tl) with statistical 
errors: (a) (top) SPI spectrum with powerlaw fit, (b) (bottom) 
ISGRI spectrum with same model as (a). The fit parameters are 
given in Sect. 14.151 The significance of the emission above 200 
keV is 3 cr. 



and Norm = 0.63+0.04. There is no evidence for a cutoff. The 
powerlaw index is still in agreement with our SPI measurement. 

4.16. GRS 1758-258 

The LMXB GRS 1758+258 is one of the well known micro- 
quasars and regarded to be similar to IE 1740.7-2942. It is also 
a black hole candidate. For a recent review see e.g. Pottschmidt 
et al. ( l2006l ). The source is bright and can show outburst which 
are however mostly too short for SPI's sensitivity. B08 and our 
analysis assume it to be a constant emitter In the low-hard state, 
which the source mostly assumes, the high energy spectrum can 
roughly be described by a hard powerlaw with an exponential 
cutoff around 100 keV. In terms of physical models, thermal 
Comptonisation models and models involving black body emis- 
sion from a disk have been applied. The spectrum above 200 ke V 
was hardly studied. 

We fit our SPI spectrum initially with a powerlaw with ex- 
ponential cutoff and obtain a reduced of 1.7 for 14 d.o.f. (Fig. 
[TTh ) and well constrained best fit parameters F = 1.54+0.07, 
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Fig. 17. The unfolded spectra of GRS 1758-258 derived from 
dataset 10 (see table[T]l with statistical errors: (a) (top) SPI spec- 
trum, and (b) (bottom) ISGRl spectrum, both fit by powerlaw 
with exponential cutoff. The fit parameters are given in Sect. 

KlE 

^cutoff = (178+28) keV, Norm = 0.12+0.03, in agreement with 
previous measurements. 

The same spectral fit applied to the ISGRI data results in 
best fit parameters F = 1.36+0.02, ^cutoff = (123+5) keV, Norm 
= 0.059+0.004, and reduced x'- = for 12 d.o.f., in marginal 
agreement with the SPI result. 

In both cases, the fit quality is not very good which may be 
the result of assuming the source to be constant (but we chose to 
follow B08 consistently). In the SPI spectrum, there is a hint of 
a hard tail in the residuals. In fact, we achieve a fit improvement 
(reduced;^'^ = 1.6 for 13 d.o.f.) if we freeze the cutoff energy and 
add a powerlaw as a second component. However, the chance 
probability of the improvement is 0.2. 

4.17. Ginga 1826-24 

The LMXB Ginga 1826-24 was discovered in 1988 and initially 
found similar in spectrum and variability to the BHC Cyg X- 
1. Not much later, the relatively low cutoff energy 58 keV) 
of its powerlaw spectrum above a few ten keV and the obser- 
vation of X-ray bursts led to the classification as a NS. The ex- 
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Fig. 18. The unfolded spectra of Ginga 1826-24 derived from 
dataset 10 (see table[T]i with statistical errors: (a) (top) SPI spec- 
trum fit with two-component model consisting of a powerlaw 
with exponential cutoff and a second powerlaw, (b) (bottom) 
ISGRI spectrum fit with powerlaw with exponential cutoff. The 
fit parameters are given in Sect. 14.171 



ceptional precision of the reoccurrence of the bursts (period ^ 
3.5 h) also led to the nickname "clocked burster". For a recent 
review of the source properties see Thompson et al. (2005). The 
time-averaged high-energy spectrum up to about 300 keV was 
measured by Strickman et al. (I1996al l and modelled as a pow- 
erlaw with exponential cutoff. They found F - 1.59+0.02 and 
£cu,oif = (57+8) keV. 

Fitting our SPI spectrum of Ginga 1826-24 with a powerlaw 
with exponential cutoff, we obtain a reduced of 1.2 for 14 
d.o.f. and best fit parameters F - 1.7+0.1, /Scutoff - (60+7) keV, 
Norm - 0.36+0.10 in good agreement with Strickman et al. 
(I1996al l. 

The residuals show evidence for a hard tail above 200 keV. 
Accommodating the possibility of such a tail by adding a pow- 
erlaw component to the model, we obtain from the fit to the 
SPI data (Fig. [Hi) a reduced of 0.5 for 12 d.o.f. and best 
fit parameters Normi = 0.24+0.09, F2 = 0.8+0.69, Norma (at 
200 keV) = (2.1 + 1.8) x lO-'^cm-^s-'keV '. The second pow- 
erlaw component becomes dominant above ca. 300 keV. 
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Fig. 19. The unfolded spectra of GRS 1915+105 derived from 
dataset 1 1 (see table[T]l with statistical errors: (a) (top) SPI spec- 
trum fit with two-component model consisting of powerlaw with 
exponential cutoff and second powerlaw, (b) (bottom) ISGRI 
spectrum fit by same model as used in (a). The fit parameters 
are given in Sect. 14.181 



Applying the same modelling to the ISGRI data (Fig. [TSb). 
we obtain for a simple powerlaw with exponential cutoff a 
reduced of 1.5 for 12 d.o.f. and best fit parameters T - 
1.78+0.04, £cutoif = (69+3) keV, Norm = 0.42+0.05. The resid- 
uals don't show any posibiUty for improvement. 

We conclude that our time averaged spectra agree well with 
previously published results. SPI and ISGRI data are compatible 
with each other The evidence for a hard tail above 200 keV is 
weak even though BOS still find 25 mCrab at 200-600 keV. We 
find less than half of that flux. 

4.18. GRS 1915+105 

The LMXB GRS 1915+105 is one of the most studied galactic 
black holes. For a review see e.g. Fender & Belloni (2004) or 
Rodriguez et al. ([2008al l. The object is a microquasar, the first 
to exhibit superluminal motion in its jets, with a relatively large 
black hole mass of ^ 14 Mq, a large binary period of 33.5 d, 
and an extremely complex time variability which is attributed 



to accretion disk instabilities (Greiner et al. 119961 Done et al. 
120041) . 

GRS 1915+105 was extensively monitored by INTEGRAL. 
A detailed study of two years of the INTEGRAL observations 
was published in Rodriguez et al. (2008a 2008b). They analyse 
ISGRI spectra up to 300 keV and distinguish several spectral 
states which they fit with a number of multi-component mod- 
els involving for the high energy part thermal Comptonisation 
(Xspec model COMPTT) with seed photon temperatures be- 
tween 0.7 and 2.1 keV and a powerlaw with index between 2.0 
and 2.9 . 

Our time-averaged spectrum of the source is a superposition 
of all spectral states. As for the previous sources we initially try 
to fit a powerlaw with exponential cutoff. This results in a high 
reduced of 2.4 for 14 d.o.f., a photon index of 3.30+0.02, and 
a cutoff energy far beyond the upper end of our energy range. 
The residuals indicate that a more sophisticated model is needed. 

Next we try a powerlaw with exponential cutoff with an addi- 
tional powerlaw. This results in a reduced of 0.6 for 12 d.o.f. 
and badly constrained best fit parameters Fi = 1.7±1.7, EcmoS 
= (13+9) keV, Fz = 2.8+0.2, Normi = 3.7+16, and Normz (at 
200 keV) = (1.4 + 0.1) x lO^^cm-^s-'keV-', i.e. a dominant 
pure powerlaw with a weak additional low-energy component 
becoming negligible at a few ten keV. 

Applying the same model to the ISGRI data results in 
F2 = 2.87+0.05 and Normj (at 200 keV) = 1.17 ± 0.04 x 
lO^^cm^^s *keV ' for the dominant component and Fi = 
0.8+1.0, ficutoif = 8+2 keV, Normi = 0.09+0.25 for the weak 
low-energy component. The reduced is 0.7 for 10 d.o.f. The 
agreement with the SPI result is good. 

4.19. CygX-1 

The HMXB Cyg X-1 consists of a black hole candidate of 
high but uncertain mass and an OB supergiant (see Gies et al. 
2003 Ziolkowski 2005 'j and references therein). It is one of 
the brightest hard X-ray sources and has been studied exten- 
sively by many authors in all wavebands. For a brief review see 
e.g. Szostek & Zdziarski (I2007I I and Shaposhnikov & Titartchuk 
(j2006|) and references therein. The radio emission shows evi- 
dence of a moderately relativistic jet. The object has therefore 
recently been called a microquasar (e.g. Russell et al. l2007l l. 

Long-term monitoring has revealed two main spectral states 
like in other BH binaries, the low-hard state and the high-soft 
state, between which the source transits. It spends the domi- 
nant fraction of time in the low-hard state (Zhang et al. 1997] 
Hjalmarsdotter et al. 2008). Our time-averaged spectrum of 
Cyg X- 1 should therefore be more similar to the low hard state 
which is roughly a hard powerlaw (F between 1.4 and 2.1) with 
a break at £ > 50 keV (Cadolle Bel et al. l2006l) . 

The X-ray spectra of individual states have been modelled up 
to 300 keV by several authors with a variety of combinations of 
disk black body and thermal Comptonisation with additional Fe 
fluorescence line and reflection components (see again Cadolle 
Bel et ai aom . 

Lacking the soft X-ray information and averaging over all 
source states, we fit our SPI spectrum of Cyg X-1 initially with 
a simple powerlaw with exponential cutoff and obtain a very 
large reduced of 5.3 for 14 d.o.f. The best fit parameters (F = 
1.65+0.01, ^cutoff = (180+5) keV, Norm = 2.3+0.1) agree wefl 
with the expectations from previous measurements and, up to 
about 300 keV, the model describes the data reasonably well but 
there is a clear and very significant hard tail above 300 keV. 
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Cyg X-1 (SPI) 
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Energy (keV) 



Cyg X-1 (ISGRI) 




Energy (keV) 

Fig. 20. The unfolded spectra of Cyg X- 1 derived from dataset 
12 (see table [T]i with statistical errors: (a) (top) SPI spectrum 
fit with two-component model consisting of a powerlaw with 
exponential cutoff and a second powerlaw, (b) (bottom) ISGRI 
spectrum fit by same model as in (a) . The fit parameters are 
given in Sect. 14.191 



This hard tail in the low-hard state spectrum of Cyg X-1 
was already observed up to about 5 MeV with COMPTEL by 
McConnell et al. (2000). They did not observe an indication of 
a cutoff. It was also already studied with INTEGRAL by Cadolle 
Bel et al. (l2006l l and Malzac et al. (|2006] |. 

Adding a powerlaw to our model to accommodate the hard 
tail, but leaving all parameters free, we obtain a much bet- 
ter fit (Fig. l20h ) with reduced 0.9 for 12 d.o.f. and resid- 
uals which don't show any sign of systematic deviations. The 
best fit parameters are Fi - 1.3+0.2, -Ecutoff = (95+15) keV, 
Normi = 0.43+0.35, F2 = 2.11+0.1, Normj (at 200 keV) = 
(9.9 ± 1.8) X 10-^^ cm-^s-ikeV-i. 

Applying the same model to the ISGRI data (Fig. l20b). we 
obtain a reduced of 1.5 for 10 d.o.f. and best fit parameters 
Fi = 1.44+0.14, £eutoff = (115±18)keV,Normi = 0.9+0.5, F2 = 
2.2+0.14, Norma (at 200keV) = (3.6+2.3)xlO-5 cm-VkeV-', 
i.e. the best fit parameters are in agreement with those from SPI 
but are not as well constrained because the ISGRI spectrum is 
lacking high energy sensitivity. 
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Fig. 21. The unfolded spectra of Cyg X-3 derived from dataset 
13 (see table [1]) with statistical errors: (a) (top) SPI spectrum 
fit with two-component model consisting of a powerlaw with 
exponential cutoff and a second powerlaw, (b) (bottom) ISGRI 
spectrum fit by same model as in (a) . The fit parameters are 
given in Sect. 14.201 

To further compare with the COMPTEL data, we fit a pow- 
erlaw to our SPI data above 500 keV. We find an index of 
1.9+0.5 which is much harder than the index 3.3+0.4 measured 
by McConnell et al. (2000) in the neighbouring energy range. 
One explanation of this difference could be source variability, 
which has indeed been suggested to be present as a variable 
broad fine feature around 1 MeV (McConnell et al. 19891 12OOOI 1. 

We conclude that even though our spectra have extremely 
high statistical significance, it is possible to model them with 
a phenomenological model with only five parameters consist- 
ing of a cutoff powerlaw and second powerlaw. The cutoff 
powerlaw can be interpreted as a contribution from thermal 
Comptonisation as already suggested by many previous authors 
(see above) while the hard tail is most likely a product of the jet. 

4.20. Cyg X-3 

Cyg X-3 is a HMXB with an unusually short binary period 
of 4.8 h and a Wolf-Rayet star as an unusual stellar compan- 
ion. Despite its brightness and the intensive study at all wave- 
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bands over four decades, recently published articles still declare 
many of its properties as uncertain, in particular the nature of the 
compact object contained in it. Szostek & Zdziarski (I2008I I and 
Hjalmarsdotter et al. (|2P08 ), H08 in the following, give a brief 
review of the literature. They also argue that the compact object 
in Cyg X-3 could be a quite massive black hole. The object has 
also been called a microquasar because of its radio jets (Marti 
et al. I2007I I. At lower energies, the intense stellar wind from 
the companion star leads to strong absorption features. Above 
25 keV, however, the effect on the spectral shape is small (HOB). 

Cyg X-3 has a binodal distribution of soft X-ray count rates. 
The distribution of its spectral hardness, however, only shows 
one peak at low hardness with a long asymmetric tail to higher 
hardness (H08). The authors note that the soft X-ray hardness 
in the case of Cyg X-3 may be a bad indicator of the state of 
the source because of its strong absorption. The source probably 
stiU has a separate hard state at higher energies which it occupies 
roughly two thirds of the time. 

The time-averaged spectra we are studying here should 
therefore be a superposition of spectra from all states, but more 
similar to the harder spectra. 

We initially fit our SPI spectrum of Cyg X-3 with a pow- 
erlaw with exponential cutoff. This results in a relatively high 
reduced of 1.8, an index F - 3.26+0.04 and a cutoff beyond 
the upper end of our energy range. The residuals suggest a hard 
tail beginning at about 100 keV. Adding a powerlaw as second 
component, we obtain a reduced of 0.8 for 12 d.o.f. (Fig. 
l2Th ) but the fit parameters are poorly constrained (Fi = 2.0+0.7, 
£cutoff=(20+10) keV, Normi = 3.5+6.6, F2 = 2.4+0.3, Normj 
(at 200 keV) = (9.0 + 1.1) x 10"^ cm-^s-ikeV-') because the 
best fit cutoff energy is outside our energy range. The powerlaw 
component without cutoff constitutes the dominant part of the 
emission above a few 10 keV. There is no evidence for a cut- 
off at higher energies, but our data points become insignificant 
above 350 keV. 

An unusually low cutoff energy implying a most likely 
plasma temperature of only 4 keV was also found in the 
INTEGRAL data analysis by H08 when assuming a thermal 
Comptonisation model. 

When inspecting our ISGRI spectrum of Cyg X-3, we find 
evidence for enlarged systematic errors which are probably 
caused by the crowded field around the object including the 
nearby and extremely bright Cyg X-1. We increase the statis- 
tical errors by 0.5 % of the flux. With this change, the fit of the 
same two-component model as for the SPI data with /Scutoff fixed 
to 20 keV results in a reduced of 1.7 for 11 d.o.f. and best 
fit parameters F, = 2.00+0.08, Normj = 3.4+0.7, F2 = 2.7+0.2, 
Normz (at 200 keV) = (5.3 + 0.5) x lO-*" cm-^g-ikeV-' (Fig. 
12Tb ) in agreement with the SPI results. 
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Systematic errors of flux measurements in crowded 
fields 
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As described in Sect. 13.31 we derive our systematic error esti- 
mate from the study of our Crab spectrum, i.e. from a source 
which is not in a crowded field as opposed to many of the X-ray 
binaries studied in this work. We need to verify to which ex- 
tent our analysis is subject to additional systematic errors when 
applied to sources in crowded fields. In order to assess this pos- 
sibility, we have randomly added a source position "Control!" 
to the input catalog for one of our most complex fields (dataset 
9) near GX 354-0 and 4U 1700-377 (see Table[I]l. This position 
constitutes the worst case for our analysis as it is close to vari- 
able bright sources in a location in the galactic plane less than 



100 200 

Energy (keV) 



> Control 1 (ISGRI) 




100 200 
Energy (keV) 



Fig. 22. The unfolded spectra for the control position 
"Controll" at RA = 259.5°, Dec = -40.35° derived from dataset 
9 (see Table[T]l with statistical errors as for all other sources: (a) 
(top) SPI spectrum with powerlaw fit, (b) (bottom) ISGRI spec- 
trum with powerlaw fit . The fit parameters are given in Sect. 



15° from the galactic center where also significant diffuse emis- 
sion is present. Furthermore, it is about 10 arcmin away from 
IRXS J171824.2-402934, an X-ray burster capable of very rare 
4 minute bursts with peak intensities up to 1.3 Crab and a ther- 
mal spectrum peaking at 5 keV detectable up to 28 ke V (Kaptein 
et al. I2007I I. As a control case, this position therefore also gives 
an estimation of the influence of variable and on average faint 
sources which are not included in our input catalog. 

The resulting SPI spectrum (Fig.l22h) does not contain any 
point with a flux significantly different from zero. All points have 
less than 2 cr significance. The only indication of any systematic 
signal which could fake the presence of a source is in the first bin 
near 30 keV where a level of 3 mCrab is reached. If we naively 
fit the points with a powerlaw, we obtain a very steep spectrum 
with best fit parameters F = 5.5, Norm - 2075. 

Repeating the same procedure for the corresponding ISGRI 
spectrum (Fig. 122b) results in best fit parameters F = 1.0, Norm 
- 1.3x10"^'. Also here all points are consistent with zero. The 
most significant point is near 300 keV and has 2 cr with a nomi- 
nal flux of approx. 80 mCrab. 
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The fluxes derived from these powerlaws can serve as ad- 
ditional estimations of systematic errors on the flux values mea- 
sured by us. They are given in Table|5]in the row labelled "Estim. 
Syst. Error" with negative sign in order to indicate that the val- 
ues need to be subtracted from the other flux measurements in 
the table. 



5. Summary and conclusions 

Tables|2][3]lll and|5]summarise our spectral measurements. They 
provide a catalog of the time-averaged spectra of the brightest 
known soft gamma-ray sources. They also give an impression of 
the variety of possible spectral shapes for different source types, 
especially for X-ray binaries (XRBs) which represent more than 
half of the sources. Furthermore, the tables permit to compare 
the response of SPI and ISGRI. In particular Tables|2]and[3]per- 
mit to study the agreement of the corresponding measurements 
in terms of spectral shape while Table |5] can be used to compare 
the flux normalisation of the two instruments in several energy 
bands. In the latter table we also give the results from BOS for 
those energy ranges where they are available. 

5.1. Agreement of SPI and ISGRI spectra 

The study of 20 bright sources with SPI and ISGRI using nearly 
identical datasets permits us to carry out a quantitative compar- 
ison of the response of the two instruments and test the quality 
of the spectra obtained with our new spimodfit-hased analysis 
method. 

The comparison of the absolute flux normalisation should be 
carried out using a steady source such as the Crab. Our Crab 
flux measurements are given in Table |4] While the flux ratio 
SPI/ISGRI is consistent with a constant value of about 1 . 12 up to 
200 keV, we observe a clear deviation at higher energies suggest- 
ing systematic differences in flux calibration of additional 18 % 
at 200 - 600 keV (see Sect. lOT i. 

If we correct for the differences in absolute flux calibration 
by normalising on the Crab flux, we obtain the flux values in 
Table |5] Plotting the SPI values versus the ISGRI values from 
this table, we obtain Figs.l23h. b, and c. In all three energy ranges 
up to 600 keV, these plots show clearly a good agreement be- 
tween the flux values from the two instruments for all sources 
with fluxes higher than a few ten mCrab. Only for the faintest 
sources and only below 100 keV do we see obvious systematic 
deviations from the linear relationship. Here the SPI fluxes tend 
to be lower than those from ISGRI. The origin of this is prob- 
ably a combination of two causes: (1) systematic differences in 
the handling of source variability by the analysis software, and 
(2) systematic differences in the background determination. 

5.2. Agreement of our SPI measurements with those by 808 

To further verify our analysis method, we compare our SPI 
fluxes with those published by BOS who give their flux values 
as fraction of their Crab flux. All values are given in Table |5] 
Plotting our values versus those from BOS (Figs. l24h . b, and c) 
we see that the flux values agree well in all energy bands for 
sources brighter than 20 mCrab. Below this flux level, our flux 
values seem to be systematically slightly lower. At least some 
of the differences can be attributed to source variability as our 
dataset is not exactly the same as that used by BOS. 

BOS defined the annihilation (positronium) emission be- 
tween 300 keV and 511 keV as an individual component in 



their source model and extracted its intensity as a result of their 
fit. Since this emission is extended (to first order a S° FWHM 
Gaussian centered on the galactic centre) on a scale similar to 
the SPI field-of-view, it is in our analysis absorbed in the time- 
variable isotropic component of our background model which 
does not have a separate annihilation component. This is con- 
firmed by the fact that our fluxes in the range 200 keV - 600 keV 
agree very well with those by BOS (see Fig.l24li. also for sources 
near the galactic centre. 

5.3. Emission above 200 keV 

We detect all of the 20 sources in our sample independently in 
the 25-100 keV and in the 100-200 keV band both with SPI and 
with ISGRI. Above 200 keV we find evidence for emission of 
more than 3 cr in S cases. In the 600-1000 keV band, we find 
such evidence only for two sources. 

XRBs constitute 70% of our source sample. Table [3] sum- 
marises our measurements of the spectra of these 14 XRBs. They 
are not a random sample but the 14 brightest known X-ray bina- 
ries at 200 keV. Four of them are HMXBs, the remaining are 
LMXBs. Ten of them contain black hole candidates (eight of 
which are classified as microquasars), the remaining four con- 
tain neutron stars. 

5.4. Conclusions 

We present spectra for the 20 brightest soft gamma-ray sources 
averaged over timescales of years and give flux values in four 
energy bands between 25 keV and 1 MeV. 

For SPI, these measurements were derived using a new tool 
(spimodift by Halloin & Strong ( I2007l l) which can coherently 
treat data from the whole SPI mission by performing a single, 
large maximum likelihood fit. 

The agreement between our SPI and ISGRI measurements 
is good if we correct for systematic calibration differences by 
normalizing on the Crab spectrum. Our SPI flux measurements 
also agree well with those by Bouchet et al. (200S) (BOS). Our 
Crab spectrum shows the also previously observed break near 
100 keV where the spectral index softens from 2.1 to 2.2. 

Up to 200 keV, all 20 sources in our sample are detected 
independently in two adjacent bands. At 200-600 keV we de- 
tect eight sources, and at 600-1000 keV we detect two sources. 
The spectra we find agree well with the results from previous 
publications. For six of the 14 XRBs in our sample (not count- 
ing the borderline case GRS 175S-25S), we find evidence for 
a hard powerlaw-component which becomes dominant above 
the cutoff energy of the thermal Comptonization component, 
i.e. above about hundred to a few hundred keV. For GX 339- 
4, GRS 1915-1-105, Cyg X-1, and Cyg X-3 such hard tails have 
already been noted by previous studies. For OAO 1657-415 and 
Ginga 1S26-24, our study provides the first weak indication of 
such emission. 
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Fig. 23. Comparison of our flux values obtained from SPI with 
those from ISGRI for the energy bands (a) (top) 25-100 keV, 
(b) (middle) 100-200 keV, (c) (bottom) 200-600 keV. In each 
band, a linear function was fit to the points and fit parameters 
are shown in the figures. The flux values are also given in Table 
|5] Only the statistical errors are shown. They take into account 
the uncertainty of the Crab flux. See Sect. 15.11 
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p1 1.02610.01122 



10 10^ lO" 

BouchetOS Flux [mCrab] 

Fig. 24. Comparison of our flux values obtained from SPI with 
those by BOS (also obtained with SPI, largely using the same 
data but with a different analysis method) for the energy bands 
(a) (top) 25-100 keV, (b) (middle) 100-200 keV, (c) (bottom) 
200-600 keV. In each band, a linear function was fit to the points 
and fit parameters are shown in the figures. The flux values are 
also given in Table|5] Only the statistical errors are shown. They 
take into account the uncertainty of the Crab flux. See Sect. [ 
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Table 2. Summary of the spectra obtained in this study (as shown in the figures of Sect.|4]i. Part la: spectral fit parameters for sources 
other than X-ray binaries. 



Object 


Instru- 
ment 


r. 


S 

[keV] 


pectral 
[keV] 


Model Fit R 

r2 


esults" 

Norm 

[keV-'cm-^s-'l 


^Vd.o.f. 


Crab 


SPI 
ISGRI 


2.11±0.01 
2.119±0.006 


- 


81 + 11 
96+4 


2.20+0.01 
2.36+0.02 


11.6+0.5 
10.7+0.3 


0.5 
1.2 


Vela Pulsar 


SPI 
ISGRI 


2.13±0.10 
2.04 + 0.04 


- 


- 


- 


0.12+0.05 
0.056 + 0.008 


0.55 
1.4 




SPI 
ISGRI 


2.10±0.17 
2.03±0.06 








0.4+0.25 
0.24+0.05 


1.4 
1.4 


NGC 4945 


SPP 
ISGRI* 


1.7±0.17 
1.56+0.05 


150 
150 






0.06+0.04 
0.033+0.006 


0.4 
0.6 


Cen A 


SPI 
ISGRI 


1.89+0.05 
1.80+0.02 








0.20+0.04 
0.117+0.009 


1.2 
0.9 


Swift J1656.3-3302 


SPI 

ISGRF 


1.5+0.2 








0.004+0.004 


1.0 



" Spectral parameters not used in the model for a particular object are marked as 
* Fit model also included the strong photoelectric absorption with A'h fixed to S.SxlO^'^cm"^. 
ISGRI spectral fit failed, see Sect.l48l 



Table 3. Summary of the spectra obtained in this study (as shown in the figures of Sect.|4|l. Part lb: spectral fit parameters for X-ray 
binaries. 



Object 


Instru- 
ment 


r, 


^cutoff 

[keV] 


S 


sectral Model Fit 
Norm 

[keV-'cm-^s-'l 


Results" 

Norm2 
[lO-'^keV-'cm-'s-'] 


^Vd.o.f. 


F-test 
prob.* 


XTE 11550-564^ 


SPI 
ISGRI 


1.36+0.09 
1.00+0.05 


206+50 
108+9 




0.031+0.010 
0.009+0.002 




1.6 
1.3 




4U 1630-47 


SPI 

ISGRI'' 


3.3+0.2 
2.3+0.2 






33+21 
0.043+0.038 




0.6 
1.0 




OAO 1657-415" 


SPI 
ISGRI 


1.23+0.07 
1.42+0.31 


23.9+3.2 
24.7+3.6 


1.27+0.58 


0.16+0.04 
0.23+0.21 


2.1 + 1.0 


0.7 
0.4 


0.015 


GX 339-4" 


SPI 
ISGRI 


1.58+0.11 
1.71+0.10 


100 
100 


0.97+0.63 
1.0 


0.081+0.026 
0.038+0.011 


3.8+2.4 
1.9+0.8 


0.8 
0.7 


0.0035 


4U 1700-377" 


SPI 
ISGRI 


2.00+0.08 
2.31+0.17 


47.5+3.9 
74.1 + 14.1 




3.14+0.71 
6.4+3.2 




1.6 
0.9 




IGR 117091-3624" 


SPI 
ISGRI 


1.75+0.20 
1.93+0.04 






0.016+0.013 
0.034+0.005 




0.4 
1.2 




GX 354-0" 


SPI 
ISGRI 


3.26+0.06 
3.53+0.02 






26.4+5.9 
43.3+3.0 




0.7 
1.0 




IE 1740.7-2942" 


SPI 
ISGRI 


1.80+0.06 
1.44+0.04 


113.4 
113.4+7.6 




0.28+0.06 
0.040+0.004 




1.9 
0.7 




IGRJI7464-3213 


SPI 
ISGRI 


2.82+0.14 
2.49+0.02 






2.1 + 1.1 
0.63+0.04 




1.1 
0.9 




GRS 1758-258" 


SPI 
ISGRI 


1.54+0.07 
1.36+0.02 


178+28 
123+5 




0.12+0.027 
0.06+0.004 




1.7 
1.9 




Ginga 1826-24" 


SPI 
ISGRI 


1.51+0.14 
1.78+0.04 


48.4+7.8 
68.5+3.5 


0.82+0.69 


0.24+0.09 
0.42+0.05 


2.1+1.8 


0.5 
1.5 


0.002 


GRS 1915+105" 


SPI 
ISGRI 


1.7+1.7 
0.3+1.1 


12.5+9.2 
7.6+2.1 


2.76+0.18 
2.87+0.05 


4+15 
0.09+0.25 


1.4+0.1 
1.17+0.04 


0.6 
0.7 


9.2x10"" 
7.8x10"'' 


CygX-1" 


SPI 
ISGRI 


1.28+0.19 
1.44+0.14 


95.9+15.5 
115+18 


2.11+0.09 
2.20+0.14 


0.43+0.35 
0.90+0.53 


99+18 
36+23 


0.9 
1.5 


9.4x10-" 


Cyg X-3" 


SPI 
ISGRI 


2.0+0.7 
2.00+0.08 


20+10 
20 


2.4+0.3 
2.6+0.2 


3.5+6.6 
3.4+0.7 


9.0+1.1 
5.3+0.5 


0.8 
1.7 


3.5x10-' 



" Spectral parameters not used in the model for a particular object are marked as 

* F-test probability that the fit improvement achieved by the inclusion of a hard tail in addition to a cutoff powerlaw was by chance, i.e. that the 
hard tail is not real. Note that since ISGRI has lower sensitivity above 500 keV and none above 700 keV, we can in most cases not test for the 
presence of the additional component but only check the compatibility with the ISGRI data. 

" Parameters of alternative fit models are given in Sect.U 

'' Spectrum very probably flawed by source contamination from IGR J16358-472. 
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Table 4. Summary of the spectra obtained in this study. Part II: Our measurements of the Crab flux with SPI and ISGRI assuming 
the spectral shape given in Table|2]for energy redistribution. 



Instrument 


25-50 keV 


Crab^ 
50-100 keV 


febula Flux [10" 
100-200 keV 


cm -s J 
200-600 keV 


600-1000 keV 


SPI 
ISGRI 


156.5+0.4 
140.2+0.4 


72.18+0.15 
64.38+0.16 


32.14+0.12 
27.46+0.09 


17.81+0.27 
13.47+0.25 


3.47+0.15 


Ratio SPI/ISGRI 


1.116+0.004 


1.121+0.004 


1.170+0.006 


1.322 + 0.032 





Table 5. Summary of the spectra obtained in this study. Part III: Source fluxes [mCrab]. 



Object Name 










Flux [mCrab]'''^ 












25-100 keV 




100-200 keV 




200-600 keV 




600-1000 keV 




SPI 


ISGRI 


SPI(B08) 


SPI 


ISGRI 


SPI(B08) 


SPI 


ISGRI 


SPI(B08) 


SPI 


Crab 


1000.0+2.0 


1000.0+2.0 1000.0+2.0 


1000.0+3.4 1000.0+3.5 


1000.0+3.7 


1000.0+11.8 


1000.0+18.6 1000.0+8.3 


1000.0+43.2 


Vela Pulsar 


10.0+0.4 


7.2+0.1 


9.8+0.7 


11.6+1.3 


10.2+0.9 


12.6+2.1 


5.6+5.3 


0.0+19.0 


13.9+4.9 


23.0+18.5 


NGC4151 


35.4±2.4 


32.7+0.7 


31.3+2.1 


43.8+7.4 


34.6+2.5 


22.9+6.8 


28.2+36.4 


0.0+132.5 


46.9+14.1 


101.4+111 


NGC 4945 


18.7+1.1 


16.4+0.3 


18.3+1.1 


21.4+3.7 


21.5+1.2 


14.3+4.0 


0.5+27.1 


0.0+53.2 


20.2+8.3 


5.8+53.8 


Cen A 


40.7+1.1 


36.0+0.3 


40.5+1.3 


59.7+3.7 


53.7+1.2 


56.9+4.6 


55.9+20.7 


87.7+24.2 


70.1+9.8 


12.9+48.7 


XTEJ1550-564 


36.1+0.6 


41.6+0.3 


23.5+0.7 


55.2+1.7 


70.2+0.9 


30.2+2.5 


48.8+8.4 


51.5+12.1 


18.8+5.5 


6.7+24.4 


4U 1630-47 


36.5+1.6 


1.8+0.2'' 


30.5+1.3 


19.7+6.0 


0.8+0.3'' 


21.7+24 


2.2+23.3 


19.9+13.2" 


21.4+5.3 


0.0+164.2 


Swift J1656.3-3302 


3.7+2.2 


d 


6.0+1.1 


8.0+3.4 


d 


5.0+1.9 


15.2+7.8 


d 


17.6+4.0 


24.6+12.9 


OAO 1657-415 


66.6+0.6 


51.5+1.3 


51.3+0.7 


15.0+1.3 


7.6+0.6 


10.3+2.2 


15.0+8.0 


0.0+24.9 


13.8+4.8 


2.3+26.6 


GX 339-4 


34.1+0.5 


11.5+0.2 


38.4+14 


32.6+1.4 


12.1+0.5 


35.6+2.4 


30.1 + 13.8 


0.0+30.3 


30.0+5.2 


66.1+274 


4U 1700-377 


171.7+0.8 


162.3+3.6 


146.3+14 


36.0+1.0 


40.9+1.6 


35.5+2.2 


13.6+8.1 


23.4+17.2 


6.4+4.8 


0.0+62.7 


IGR J 1709 1-3624 


5.2+0.8 


6.1+0.1 


8.9+0.6 


7.2+5.0 


6.3+0.3 


5.6+2.0 


18.5+8.8 


2.7+13.4 


9.8+44 


8.7+19.6 


GX 354-0 


34.2+0.5 


23.4+0.1 


24.4+0.5 


10.5+1.1 


5.4+0.2 


7.6+1.8 


0.7+10.0 


1.0+10.5 


7.9+4.1 


7.4+214 


IE 1740.7-2942 


52.8+0.9 


32.7+0.1 


51.1 + 1.2 


43.0+2.3 


35.0+0.3 


48.3+1.6 


18.5+10.9 


22.1+6.1 


39.2+3.6 


7.8+34.9 


IGR J 17464-32 13 


13.3+0.6 


15.1+0.1 


15.6+0.7 


6.9+1.1 


9.9+0.2 


11.2+1.6 


2.2+4.4 


6.8+6.3 


9.9+3.6 


41.6+23.1 


GRS 1758-258 


68.8+0.7 


66.6+0.2 


70.8+0.9 


84.6+2.1 


82.2+0.5 


83.4+1.6 


41.1+7.2 


49.3+6.0 


47.4+3.6 


51.2+18.0 


Ginga 1826-24 


75.8+0.5 


73.1+0.3 


74.6+0.5 


30.3+0.9 


34.2+0.4 


33.6+1.9 


18.4+5.7 


0.0+19.8 


24.5+3.9 


38.3+224 


GRS 1915+105 


189.1+0.8 


187.8+0.5 


150.1+0.9 


73.2+2.2 


60.5+0.5 


61.8+2.3 


38.8+8.7 


40.9+10.2 


38.6+5.2 


53.9+27.2 


Cyg X-1 


841.0+1.5 


871.6+1.7 


843.8+1.0 


873.7+3.1 


907.6+3.1 


769.9+3.1 


575.5+9.8 


581.3+16.1 


500.6+6.9 


386.4+37.3 


Cyg X-3 


117.7+0.6 


111.7+0.6 


101.2+0.8 


34.9+1.4 


25.3+0.5 


25.0+2.9 


26.4+7.6 


21.0+21.8 


14.8+6.2 


52.6+31.0 


Estim. Syst. Error* 


-0.9 


-0.1 




-0.1 


-0.4 




-0.1 


-1.2 




-0.1 



" For the columns "SPI" and "ISGRI" fluxes are expressed as fractions of the Crab fluxes for the individual instrument given in Tabled] In the columns "SPI 
(B08)" we give the fluxes obtained by Bouchet et al. 120081 (B08) for the available energy ranges (see B08 for their definition of the unit "mCrab"). 

* The statistical eiTors shown in this table do not take into account the uncertainty of the Crab flux on which the flux values are normalised. 

'' Spectrum very probably flawed by source contamination from IGR J16358-472; therefore not included in Fig. 1231 
ISGRI spectral fit failed, see Sect. 14. 81 therefore not included in Fig. 1231 

' See Sect. 1121] 



